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FIRST RESULTS OF A STUDY ON TURBULENT BOUNDARY LAYERS IN 
OSCILLATING FLOW WITH A MEAN ADVERSE PRESSURE GRADIENT 

by R. Houdeville and J . Cousteix 

Summary 

The development of a turbulent unsteady boundary layer with a /l » 
mean pressure gradient, strong enough to induce separation, is studied 
in order to complete and extend the results obtained for the flat plate 
conf igura tion . 

The longitudinal component of the velocity is measured using 
constant temperature hot wire anemometer. The investigation of the 
region where negative velocities exist is performed with a laser 
Doppler velocimeter system with BRAGG cells. 

As in the flat plate configuration, the boundary layer responds 
by a forced pulsation to the perturbation of the potential flew. Due 
to the experimental conditions, the unsteady effects are very important; 
for example, the average location of the zero skin friction point moves 
periodically at the perturbation frequency. 

The study of the longitudinal component of turbulence deals with 
its standard deviation as well as its moments of order 3 and 4 which 
give some information about the development of turbulence. 

The comparison of average velocity profiles, obtained at different 
instants in the cycle, with profiles deduced from similarity solutions 
shows a good agreement up to values of the shape parameter close to the 
separation value. Moreover, the existence of a logarithmic region 
enables a simple calculation of the maximum phase shift of the velocity 
in the boundary layer. 

To conclude, an attempt of calculation by an integral method of 
the boundary layer development is presented, up to the point where 
reverse flow starts appearing. 

1 - I NTRODUCT ION 


Numbers in the margin indicate pagination in the foreign text. 


1 


The study of fist plans boundary layars subjactad to parlodic 
oscillating flows is mainly conducted in the reduced frequency range 
which is not too high. If such experiments allow for the definition 
of a certain number of important theoretical data, they must neverthe- 
less be quite far from practical cases encountered in unsteady aero- 
dynamic problems encountered, for example, the problem of the dynamic 
jitter of a helicopter blade. 

In order to approach such a situation and also to provide experi- 
mental data required for the development of calculation methods, we 
shall study the development of a pulsated boundary layer with a positive 
mean pressure gradient which may lead to separation. The experimental 
study deals mainly with the longitudinal component of velocity as well 
as with its fluctuating components characterized by its standard devi- 
ation and its higher moments. 

2 - EXPERIMENTAL CONDITIONS 

2.1. - Experimental Device 

A turbulent boundary layer subjected to a pulsating external flow 
is developed on the floor of a test section which is rectangular: 160 x 

220 mm. Abrasive paper, glued 100 mm ahead of the inlet of the test 
section along a length of 105 mm in the collector of the wind tunnel 
allows the transition to be triggered and measured. 
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Upstream from tha study ssction is located a swirl vane, the rota- 
tion of which induces the flow pulsation by changing the total pressure 
loss of the aerodynamic circuit. 

A streamlined body, composed of rigid parts interconnected among 
each other by flexible connections, induces on the floor of the test 
Motion a longitudinal mean pressure gradient which is adjustable. 

Suction of part of the higher rate of flow prevents blocking from taking 
place toward the trailing edge of the main body, as this would induce /2 
considerable overvelocities on the floor at the test section inlet. A 
correct adjustment of the upper suction area equalizes the pressures 
of the streamlined body at the test section inlet. 

To prevent a separation of the lateral boundary limits and, hence, 
spurious tridimensional effects, a suction of these lateral boundary 
layers is produced in a region located upstream from the one where the 
unsteady separation occurs, which is the objective of this study, which 
permits the passage of the light beams from the laser anemometer. 
Visualizations of the parietal streamlines, which are performed in an 
unsteady state with a mixture of oil and black smoke have made it pos- 
sible to state that it was possible to obtain an average two dimensional 
flow over about half of the space across the test section and over a 
length ekceding 600 mm. 

As previous studies have demonstrated the detrimental effect of 
the presence of static pressure intakes on the development of a bound jury 
layer in unsteady state, none have been placed on the floor in the 
measuring region. 

Considering the extreme vibratory environment, particular attention 
has been brought to the attachment of the test section. For example, the 
translation mechanism of the measuring centers is not mechanically con- 
nected to the wind tunnel, but is attached to a stable support fixed 
directly on the ground. With these precautions taken into account, 
there is no visible evidence of vibration of the measuring sensors, nor 
of the floor of the test section from the sighting telescope. 

2.2. - Methods of measuring. Data proc e ssin g. 

The measurement of instantaneous velocities is performed in a 
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Figure 2. System of acquisition 
and treatment of data 
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d) synchronization 

e) system of data acquisition 
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g) velocity 

h) overall average 

i) time average 

j) mean square fluctuation 


standard manner with the use 
of DISA 55 M type hot wire 
anemoneter with a constant temp- 
erature. The range of application 
of the hot wire was limited to 
positive speed flows, the use 
of the laser anemometer has 
made it possible to measure the 
velocity in regions where there 
are return flows. 

The velocimeter used, of the 
Doppler type, was developed at 
ONERA [1-2]; we shall review 
briefly its operating principle. 

A coherent light beam (X * 51^.5 
mm) of 2.8 watts emitted by a 
nitrogen laser source is divided 
into two parallel beams of equal 
intensity symmetrical in relation 
to the measuring volume; this is 
carried out with the use of two 
BRAGG cells, located behind the 
beam divider, which separates 
the light ray frequencies of 200 
and 204 MHz respectively. The 
interband, about 7.85 ym, was 
determined by the measuring vol- 
ume on a photomultiplier. There 
is no special problem concerning 


space requirement in the case of our mounting, diffusion takes place 
toward the front. The processing of the Doppler signal is done using a 
DISA 55L90 computer, directly connected to a Hewlett Packard 21 MX 
computer through its digital output. The two windows having dimensions 
of 200 x 100 mm, are made of antiglare treated glass with an evenness 
of 10 X along their longest length. 



Technique of measurement analysis 
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As the flow Is periodic and turbulent, a statistical analysis of 
data is required to separate the average component of velocity from the 
turbulent part. To accomplish this each cycle is considered as the 
occurence of the same phenomenon; the averages , i . e . overall averages, are 
computed by always placing the various successive periods in the same 
phase. The phase data is provided by the rotation of the swirl vane 
which induces the periodic flows . This rotation generates, through a photo- 
electric cell, a series of impulses which always correspond to the 
specified values of the phase angle. The signal from the hot wire 
anemometer is digitized at each impulse; we have, then, the velocity 
of 24 moments in each cycle. 

In the case where the laser velocimeter is used, the procedure is 
a little more delicate, as in principle the velocity data is sampled 
and supplied randomly in time. The measured speed is the last value 
at the counter prior to the corresponding impulse. In order to effect- 
ively obtain the velocity at a given phase and not an average over 
l/24th of a cycle, which would systematically result in substantial 
errors, for example, in the rate of turbulence, the flow is seeded with 
sensor particles in order to insure very high measuring rate at the 
counter which would exceed 50,000 measurements per second, corresponding 
to an average lag between two measurements of less than l/lOOO*-* 1 of a 
cycle. In the relatively improbable case where no particle passes into 
the measuring volume between two impulses, the value computed by the 
counter is not taken into account. 

By definition, the random fluctuation of the velocity u' has an 
overall average of zero, u' is obtained by the difference between the 
overall average of the velocity U and its instantaneous value, u' may 
then be analyzed by using its standard deviation or its 

moments of higher orders. 

In order to obtain a certain accuracy for the values of U or for 

, averages must be calculated from a sufficient number of results; 
we have selected to work with 600 or 1200 cycles depending on the rate 
of turbulence. We shall estimate the accuracy thus obtained from the 
statistical theory of errors. The law of large numbers indicates in 
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fact that if X^, X 2 '>.«,X n are the values taken by a random variable X 
during n independent computations, whereas when the variable $ *</«* 

leads to a random variable which follows a standard deviation law r/*r 
where o* is the standard deviation of X; m being the mean value of 
X: r. . we have then: 



The accuracy obtained for the mean velocity U by U • £ u */" can be cal 
culatedif it is given that the probability of error is below h; with £ * 0.9, 
we have then: -iiL= 4 ,(* 

T 

Let us assume the following conditions: 

. u t o - *oo (region quite far from the separation) (1) 

h • 7 */•• 

iEL -so* n - 1200 (region near the separation) (2) 

u h - 2,4 * 

The determination of the accuracy on the turbulence rate is 
obtained in a similar manner, X^ represents now the standard deviation 
of the velocity. Let us compute the standard deviation of X^ : 

<u'V _ 


If we assume a Gaussian distribution for the turbulence, we have 


then: 


<*'*> z 3 


and 




Always with a probability of 0.9, the error made on the turbulence 
rate is given by: 


s 


It may be observed that it is independent from the turbulence level. 


For n ■ 600 and 1200, we have then h - 10% and 7% respectively. 

It would have been possible, of course, to use more points in 
order to insure a better accuracy. However, it should be observed 
that it increases as fT , or 1200 points which represent data acquisition 
during one minute for a simple point of the boundary layer. 

The computations made above for estimating the measurement accuracy 
depend on the independence of the various results; if account is taken of 
the relatively large difference between the oscillating frequency and 
the characteristic turbulence frequencies, it may be considered that 
this last condition has been accomplished, especially since the averages 
are computed every 2 cycles. 

2.3. - Laser veloclmeter - hot wire comparison . 

Before beginning the presentation of the experimental results, we 
shall make a comparison between the measurements made with the laser 
veloclmeter and the hot wire in the same conditions used for our ex- 
periment; for example, we use the BRAGG cells, even though they are not 
indlspensible in this particular case. To accomplish this we shall 
place a station at 435 mm from the test section inlet where return 
flows never occur. 

Treatment of the data processed by the laser veloclmeter 1.p the 
same as for measurements made with the hot wire; overall averages are 
computed by selecting a given phase in the cycle; as the rate of par- 
ticle flow is high for such sampling no correction is required to 
compute the averages. On the other hand only the data located in the 
standard deviation interval i 5 are taken into account; it is important 
not to reduce this interval in order to avoid the likelihood of making 
a large error in the calculation of the moments of orders 3 and 4 of 
the speed fluctuation. Finally, it is necessary to define the size 
of the measuring volume, limited by the diaphragm located in front of 
the photomultiplier, which is 0.3 mm; as the boundary layer thickness 
is more than 30 mm, the error in calculating the moments with even 
degrees by the average speed gradient in the measuring values is small /H 
compared to other causes of error, it is thus not necessary to make a 
correction in the calculation of these moments. 
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In figure 3 we shall compare the evolution as a function of time 
of the mean velocity relative to Its continuous component; this Is 
performed In the boundary layer at various distances from the skin and 
1200 cycles are used for this measurement. 



The deviations observed are small 
and are not significant; in fact, for 
practical reasons, It was not pos- 
sible to make estimates during a 
single test by placing the hot wire 
Just behind the measuring volume of 
the velocimeter. 

The turbulence intensity pro- 
files, obtained at three different 
Instants in the cycle (Fig. 4) show 
a satisfactory agreement. It should 
nevertheless be observed that the 
laser velocimeter always indicates 
an excessive turbulence rate out- 
side the boundary layer; this comes 
from a limitation in the counter 
efficiency, with the accuracy of each 
individual measurement being about 
1.5*. 


Figure 3. Anemometer comparison 

of laser-hot wire. Mean velocity. To extend our comparison further 
we present In Fig. 5 the evolution of the dissymetrical and flattening 
factors in the boundary layer at a given instant in the cycle. The good 
agreement observed would not have been obtained if we had to eliminate 
in the case of the laser velocimeter the velocity values outside of 
the interval t 2.5 or 3 r ; this was made possible due to the good 
quality of the Doppler signals and to the effectiveness of the circuits 
eliminating incorrect counter measurements. 


As the coherence of the two measuring methods used, laser velo- /£ 
cimeter and hot wire, was established within the framework of our 
study, we shall not define in the presentation which follows the ex- 
perimental results the measuring method used. Let us add only that the 


probes of the boundary limits performed up to the location x ■ 390mm 
have been measured with the hot wire, the laser veloclmeter then took 

its place. The external velocity was especially measured with the hot 

wire. 


Key: a) Laser (with Bragg 
cells) 

b) Hot wire 



laser-hot wire. Turbulence. 
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Figure 5. Anenoeeter comparison of laser- 
hot wire. Dissysetry factor ami flattening. 


3 - EXPERIMENTAL RESULTS 


3.1. - External velocity 

As we have seen the pulsation of the flow is created by the rotation 
of the swirl vane located upstream from the test section. It may be 
observed that there is a resonance frequency for which the pulsation is 
practically sinusoidal. The entire study presented is made of this 
special frequency equal to 38 Hz. 


The harmonic analysis of the external velocity limited to the 
fundamental frequency is presented in figure 6. The harmonics of 
higher order remain low, below 10% of the fundamental value, except 
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upstream from the region with the return flow (X > 550 nm) where they 
become quite high. 
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Figure 6. Harmonic analysis of 
exterior velocity. 


The continuous component of ve- 
locity, close to 27 m/s at the test 
section inlet, increases up to 30 m/s 
in the region of the constant section 
under the streamlined body; this ve- 
locity increase comes simply from the 
decrease in this section due to the 
main body. It may be observed that 
it would be considerably higher with- 
out the suction of one part of the 
flow rate passing above the obstacle. 
After X - 200 mm, the divergence of 
the channel creates a constant de- 
crease in the continuous component. 

The variation amplitude, approx- 
imately 20% of U Q at the test sec- 
tion inlet, decreases considerably as 
it moves farther upstream; it should 
nevertheless be observed that the rel- 
ative fluctuation amplitude, ie the 


one related to the local continuous component, remains essentially con- 
stant and ranges between 15 and 20%. 


The phase angle remains low with a slight increase in the decelerated 
part of the flow in addition to the Reynolds number which remains an 
essential parameter to characterize the flow, other parameters must be 
considered in the unsteady flow. Among the most important ones, we shall 
draw to attention the Strouhal number Xu> ■ wX/U m 3 the relative fluctu- 
ation amplitude of the external velocity: AU e /U e , the frequency of the 
periodic flow w and a characteristic turbulence frequency f e ■ u'/g. 

As the study of the boundary layers is carried out on the floor of the 
test section, the length X is counted from a fictive origin of the 
turbulent boundary layer determined by using a flat plate law for the 
mean flow. The table below defines the values of the various parameters. 
For the calculation of f t the characteristic velocity u* of the turbulence 
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under consideration is taken as equal to the mean value of && , 
next to the skin. This frequency is considerably higher then the 
periodic frequency of the flow, except in the vicinity of the region 
where the negative velocities exist; we may then expect a certain 
interaction of the pulsation on the development of turbulence and 
therefore more accentuated unsteady effects in this region. 

3.2. - Integral values 


Observation of the development of 
displacement thickness and of the shape 
parameter as a function of time at var- 
ious locations depending on X permits 
a first global description of the de- 
velopment of the boundary layer and of 
the beginning of separation. Figure 7 shows a substantial Increase in 
the mean displacement thickness as the location of the probe recedes; 
between the extreme locations studied, the mean value of ^ goes from 1.5 
to 18 mm. In the upstream parts of the zone under study 6^ varies in a 
steady manner as a function of time, with only a very small variation 
f *om a sinusoid; the same is the case for the shape parameter H (Fig.8). 
Furthermore, 4 and H are practically in phase opposition In relation to the 
external velocity; the boundary layer thickness passes through a maximum 
when the velocity is low. This has already been observed during studies 
performed on the flat plate configuration without the mean pressure 
gradient [3-4-9]. 

The unsteady effects seem to increase when one moves upstream, it 
may be observed in fact that the variations of H and of 6^ are 
increasingly larger. This may be explained by observing, as we shall see 
later on during the study of similarity solutions, that the pressure 
gradient parameter due to the flow fluctuation has the shape ;/6 

this term is larger as the boundary layer is thicker and the parietal 
friction is lower, i.e. when we approach the region of separation. 

It is not easy to locate at each instant the zero friction point; 
for example if we consider that c^ ■ 0 corresponding to a value of the 
shape parameter close to 2.6 (5), it may be observed that a further lo- 
cation of the zero friction point is close to 450 mm. The displacement 
amplitude of this point is quite high as at the last location of this 
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measurement (X * 604 mm), there are 
stili instants in the cycle where H 
remains less than 2.6. The presence 
of considerable harmonics may also 
be observed in this region, as well 
as for 6^ and for H; these har- 
monics disappear toward the upstream 
region* 

3.3. - Harmonic analysis of velocity 
in the boundary layer 

The study of unsteady effects on 
the boundary layer profiles may be car- 
ried out on the basis of the harmonic 
analysis of velocity. As the latter 
is almost sinusoidal, we are interested 
only in the first harmonic. 

The continuous velocity component profiles shown in figure 9 for the 
various locations depending on X do not demonstrate any special charac- 
teristics showing a large difference compared to the unsteady boundary 
layer profiles. This confirms a fact already observed in the case of 
a flat plate configuration, mainly the low effect of the nonlinear 
equation terms on the development of the average boundary layer component. 

The amplitude fluctuation profiles reduced by the external velocity 
fluctuation, are shown in figure 10 at various abscissas. The locations 
X * 100 and 196 mm correspond to a region where the mean velocity 
gradient is zero; the excessive amplitude in the boundary layer, less 
than 10*, is similar to what has been observed by various authors for 
such a configuration [3*9]. The positive mean pressure gradient begins 
to act after X * 220 mm which is shown on the amplitude profiles by con- 
siderable increase in the excess of AU e near the skin which could reach 
nearly 60% at X = 390 mm. Also observed toward the outside of the boun- 
dary layer are values of the relative amplitude considerably lower 
than 1. 

The phase angle profiles develop more as one moves upstream, i.e. 



Figure 7. Thickness of 
displacement . 
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when the mean pressure gradient is higher. For sensor locations farther 
back no more positive phase angles are observed near the skin. 



Figure 8. Shape parameter. 


3.4. - Mean velocity profiles, /g 
Study of turbulence . 

After describing the over- 
all evolution of the velocity 
profiles according to X , we shall 
now present the same profiles 
but fixed at various instants 
in the cycle, and that for two 
measuring locations. These 
profiles, as well as the tur- 
bulence intensity profiles, have 
been estimated as overall 
averages from 600 and 1200 cycles 
for the stations 240 and 604 
nun respectively (Figs. 12 and 
13). 

At station x = 24Q mm, i.e. 
at the location where the mean 
pressure gradient begins to move, 
the velocity profile varies 
moderately over the cycle, as 
the boundary layer is subjected 
to an instantaneous variable 
pressure gradient. The effect 
of this gradient is easier to 


observe on the evolution of turbulence profiles (Fig. 12), their shape 
presents at instant 21 a pattern similar to the one observed for an 
unsteady boundary layer which is slightly decelerated. 


The location of station 604 mm is very different (Fig. 13) at 
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Harmonic 
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certain instants, 17-21, the mean velocity direction is reversed 
near the skin. The return flow region, which takes up 1358 of the 
boundary layer thickness, is perfectly visible on the correjponding 
profiles. At instants 5 and 9 the pattern of the profile is character- 
istic of a retarded boundary layer. At instant 1 it assumes an inter- 
mediary shape having a less standard pattern. It may also be observed 
that the variation is smaller for the thickness of the boundary layer 
during the cycle, which seems quite surprising if one recalls that the 
displacement thickness varies almost by a factor 2. 
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Figure 12. Profiles of velocity 
and turbulence at X = 240 mm. 


The question arises here for 
this station whether or not at in- 
stant 17-21, the boundary layer is 
separated. Contrary to the 
steady case, in fact the sep- 
aration does not correspond to the 
point where the friction is reduced 
to zero and where it is very possible 
to have a negative mean velocity near 
the skin without the boundary layer 
being separated (14-5). Theoretical 
studies on laminar flows (13) or 
numerical studies for turbulent cases 
(11-5) demonstrate that the unsteady 
separation must be combined with a 
fast thickening of the boundary 
layer. Such a thickening does not 
seem to exist in the present study. 

Perhaps a visualization of the bound- 
ary layer with the use of smoke would 

make it possible to resolve the 
indeterminacy about the nature of 
the flow. 


The intensity profiles of turbulence at this station show a quite 
high maximum ranging between 11 and 1758 of U e , which is quite charac- 
teristic of an intense positive pressure gradient. It should be pointed 
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out that the high value of the turbulence rate in the external flow 
comes uniquely from the limitation of the accuracy of the laser 
anemometer. In fact, the hot wire shows an external turbulence close 
to 1% of U g at this station. Finally it may be stated that a sub- 
stantial harmonic rate on the external velocity exists and this perhaps 
is explained by a strong coupling between the potential flow and the 
boundary layer in this region. 



The study of moments of order 3 and 4 of the fluctuating component 
' of velocity contributes additional information relating to the develop- 
ment of turbulence. Figures 14 and 15 show the evolution of the /9 

dissymeterical factors S« <**'>/ < and flattening factors 
F s <*•>/< at two locations of the probe and at various instants 
in the cycle. The evolution of these two parameters makes it possible 
to define the state of the distribution function u* related to a Gauss 
curve for which S * 0 and F = 3. Near the skin, at the station 

^ m 1Q q x is mo( jerately negative 

| U but remains near 0 whereas F is about 

m ^ i ,* SI. 3. The distribution function of u’ 

- 1 must therefore not vary very much 

*1 from a Gaussian distribution. This 

— — i# has already been observed in the case 

of steady boundary layers and unsteady 
■ y y flat plate configurations (4). This 

/ Jr 1 is not the case at station X ■ 604 mm 

V i and especially at Instant 21 which 

corresponds to a value close to 4.4 

U« of the shape parameter, therefore to 

A i a re turn flow profile. Compared to 

• /y a ® auss curve » the distribution u* 

/•/-''' • near the skin is longer with a peak 

(jj'sr \\\ on the s ^ e °f negative values of the 

f velocity fluctuation. 

On the outside of the boundary 

mm 

• L_ ■ ' layer it may be observed that the 

• a 4o 

dissymetrical factor becomes strongly 



mm 


Figure 13. Profiles of velocity 
and turbulence at X * 604 mm. 


negative and that the flattening factor increases, and this is shown 
on the u' histogram by a larger peak on the side of the positive 
velocity fluctuations. This contributes some knowledge concerning the 
motion of the three zones of the boundary layer. The location of this 
zone has a random size, of which the distribution function is called 


the intermit tance factor; provided 
intermittance factor y is computed 


i 



Figure 14. Factor of 

dissymetry and flattening 

X = 100 mm. 


that some assumptions are made, the 
simply from the flattening factor 
by the relationship y * 3/FJ at a 
given instant y(y) represents the 
probability for which the zone of the 
boundary layer is farther from the 
wall. The curve variations of y 
corresponds to a mean motion of the 
boundary layer zone induced by the 
flow pulsation. The amplitude of 
this forced oscillation of the free 
zone, related to the mean thickness 
of the boundary layer, is greater 
than X = 100 mm than at the station 
where negative velocities are present. 
In the external flow, a Gaussian 
distribution of turbulence is 
found again. 

In the region where a return 
flow exists at certain instants, it is 
important to know the probability so 
that at a specified instant the 
velocity is negative. The curves on 
Fig. 16 show this probability as a 
function of the distance from the 
skin at various points of 
the maximum of the shape parameter H. 
When H increases, the zero friction 
point goes back upstream, whereas it 
displaces in the direction of the 
mean flow when H decreases. Quali- 
tatively the curves corresponding 
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to the instants where the point C f - o reverses (t ■ 21,23,25) are 
very similar to the ones obtained in the case of unsteady separation 
[6], This in no longer the case when the point C f * 0 goes back up- 
stream. It is observed that the negative velocities seem to be located 
closer to the skin at instants 13 and 15. This is more clearly obvious 
on Pig. 17 where there are plottings of the isoprobability lines for a 
negative velocity as a function of X and y at various instants. The 
negative velocities take up a thin region near the skin when the zero 
friction point goes back up t h e flow (instants 15 and 17). On the 
other hand, this region Is thicker when the point C^, ■ 0 recedes 


(instants 21-23). The slope of the 
higher. 



isoprobability lines is then much 

According to the first results 
of our study, it doesn’t seem that 
the unsteady characteristic of the — 
flow profoundly modifies the be- 
havior of the boundary layer sub- 
jected to a positive mean pressure 
gradient which is quite intense. 

For example the mean velocity pro- 
files have standard patterns, as 
is also the case for the instanta- 
neous profiles of turbulence. In 
order to define this first conclusion 
which is somewhat subjective, we 
shall now attempt to analyze the 
experimental results already shown 
through the use of similarity 
solutions . 

i4 - ANALYSIS OF EXPERIMENTAL RESULTS. 

SIMILARITY SOLUTIONS. 


Figure 15. Factor of dissymmetry '•- 1 - ' Petermlntlng the similarity 
and flattening at X * 60^4 mm. solutions 
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Figure 16. Profiles of’ probability of having 
negative velocities at X * 604 non. 
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Figure 17. Charts of isoprobability of having 
negative velocities 









The theory of connected asymptotic solutions shows that If the 
Reynolds number moves toward infinity, it is possible to distinguish two 
regions in the boundary layer: 

- The region of the skin in which the phenomena depend Inainly on the 

friction of the skin and on the viscosity and where the variables used 

are: ' , __ 

U s vA» c it j with Vf * 

- External region in which the friction is due uniquely to the turbulence 
and for which the flow is considerable for external conditions and 
especially for the pressure gradient, in this region the variables 

ere used: 

F' s itniL and 7 • -2- 
'l 

As the Strouhal numbers for flows which we consider not to be too /n 
high, we shall simply assume that the law of this skin defined in its 

state remains valid. This leads then to the overlapping regions, 
i.e. where y + moves toward infinity and n moves toward zero, at a ve- 
locity profile represented by a logarithmic law: 

u\ Ufwith X "°’ 4 ‘ (1> 

To determine the velocity profile shape in the external region we 
are making the assumption that the adverse velocity profile F' ■ Is 

a function uniquely of ,this is a similarity assumption, 6 being the 
standard thickness of tne boundary layer function of x and t: 




( 2 ) 


Since the Reynolds number moves toward infinity, the viscous friction 
can jjg omitted in the presence of the Reynolds stress - c <“'•"> 
this stress is expressed simply by using a combined length model (10) : 


with 

20 



( 3 ) 


By using the similarity assumption (2) in the turbulence 
diagram, we can write the local equation for the quantity of motion 
Mr the variable F' . She assumption of an infinite Reynolds number 
makes it possible to expand this equation in. relation to the small parameter 
Y*fE77T . We shall select only the terms of order *5* ; we thus 


have: 


V JL. 


(1) 


«tth e-- 

y*( < V *) - 

with the boundary conditions: 
to the relationship: 


rv» 

t ■ 0 and F* 


0 in n ■ 1 leading 


with 


-i-2L * -i-ii.- 

/V;, 


- 


(5) 


and equation (4) is finally put in the form: 

fc)‘ F ' L -- *- -jr- * f 1 f ' <6> 

with F * 


Ihe local equation of the quantity of motion has thus become 
a common differential equation for the adverse velocity F^ . Its 
solution does not depend only on the parameter 0* as F^ is determined 
when the solution is known. Ihe only modification relating to the 
steady case is the substitution of the parameter 20 by 0^ + 20 x . 
This means that the family of profiles obtained In the unsteady state 
is the same as that obtained for the steady case parametered by 
•using 0* + 0., and not 20.. 

v X X 
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4.2. - Analysis of the Similarity Solutions 


The parameter $* does not allow analysis of the family 
of profiles, as It Is directly related to the local gradient 
of external velocity which Is usually Incompatible with the existence 
of a similarity law which implies that 6* must be constant, therefore 
independent of x and t; this Is a condition which Is not rigorously 
fulfilled except for the boundary layer In equilibrium. In general 
It Is preferred to use as a parameter of the family of profiles a quantity 
which characterises the shape of the profiles themselves. We have 
selected the Clauser factor defined by: 

6 ■ 


By definition P' ■ (Ue - U)/Ut will have the following 
relationships: 


(7) 

G-j (M*4 )/m* 
it/S , k pj 

( 8 ) 

where P^ Is a function of 0 determined by the similarity solutions 
which 'are represented analytically by the formula: 

Ft -Wie -{V ♦ wwf*/* - 

In order to be able to compare the experimentally measured 
'profiles to the profiles deduced from the similarity solutions. It 
Is necessary to use a law of friction to determine the parameter 

t, iUTT . . This law is obtained simply by over- 
lapping the skin law for • and the law of adverse velocity 
for n ♦ 0 . The similarity solutions show that when n -*• 0 adverse 
velocity profiles follow a logarithmic law of the shaoe: 
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% l 


( 9 ) 


F's 1’i-u - . 

u* % » 


where D is a function of G, determined by these solutions. 

The . relationship (8) and the law of the skin (1) lead to the 
friction law: 

i: * Ut. » 


. ± z ± L m idii. * U m f t 

assuming (10) i * v * 


( 10 ) 


D* ■ 5*25 ♦ D - (1/x) L r Pj Is a function of 0 which may be 
computed from the similarity solutions and represented analytically 
by the expression: 

y* t 2 6 - <,« 6 * * 


(ID 

. The relationship (7) combined with the friction law (10) permits 
the determination of the Clauser parameter G of a given experimental 
profile if H and 6} are known. In figure 18 a comparison Is provided of 
the experimental profiles obtained at various Instants In the cycle 
with profiles deduced from the similarity solutions and corresponding 
to the same value of the parameter G. It may be stated that a good 
overall agreement exists, Including one for the profiles of which the 
shape parameter Is quite high, of order 2. We may then conclude 
that the unsteady profiles may be represented by a profile family 
for a parameter, even for boundary limits subjected to intense 
positive pressure gradients. 


The presence of a logarithmic region for velocity profiles /12 
Is shown in Pig. 19. 3he coefficient of friction present in the 
logarithmic law (1) was determined to have the best possible agreement. 
It may be stated that as with a steady flow, the extent of the logarithmic 


23 












region decreases as the shape parameter Increases, while the 
experimental Reynolds number remains finite. 


* ■ 


4.3. - Development of the logarithmic law. 
Calculation of the phase angle . 


The existence of a friction law (10) and of a logarithmic 
velocity profile In the overlapping region (1) permits the cal- 
culation of the phase angle for velocity In this region. For 
that we make the assumption of small perturbations, l.e. we assume 
that the amplitude of the various values present In these relation- 
ships remains small, which must be the case If the amplitude of the 
external velocity fluctuations Is not too high. We have t>*n: 

ut • ♦ 'Jktt 1 '** 

U % * Vt i. 

where only la complex. Adding these developments togethe &ith ai 
those of y» H, 0, 6 1 and 0 In the relations’ (1) (10), we obtain 
the following expression for velocity in the logarithmic region: 


v? Y« M,, (ia. 1 »< * s,*i '! 

% v 

♦ * M-i}- M >"' 


( 12 ) 


We see then that the phase angle of velocity In this region 
Is Independent of y and Is equal to the phase of y./ y. ♦ LL./tJ 

• • CJl VO | 

which Is also the phase of ^ . The velocity phase In the 

overlapping region Is therefore constant and equal to that of the 
parietal friction. 


Based on the law of friction (10), on relation (7) 
and on making an estimate of the small perturbations, we obtain 
for the phase velocity the expression: • 






( 13 ) 
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where AU®, A61, A6 representing the fluctuation amplitudes of 
U e » $1 and 6 and where the coefficients C and D do not only depend 
on the mean flow: 

C z J±. j>« 3 

X A A 

tilth 

A; (l - *-*)/** -1 

6* (4 - 

The phase angle which is thus calculated is represented by the 
arrows on Fig. 11 for several experimental profiles. It corresponds 
to a maximum of the phase angle located In the logarithmic 
region. The values of y + which correspond to this maximum are 
included, for the various profiles represented, between 30 and 100. 
Experimentally, it ir-.y be stated that the phase angle decreases near 
the skin which does not seem very concordant with the fact that the 
parietal friction phase must be equal to that of the velocity in the 
logarithmic region. Tb eliminate this difficulty we must assume that 
the flow of velocity in the sub-layer is not unique but is a function 
of time. We must in any case insist on the difficulty arising when 
measuring a phase angle requiring an accuracy of several degrees. 
Moreover, the decline in the phase occurs very close to the skin, 
in a region where it is not excluded that a certain interaction of 
the skin on the thermic field of the hot wire exists. The measurements 
computed with the laser anemometer nevertheless confirm the results 
observed. This difficulty shows the importance of a detailed 
study of the sub-layer, unfortunately impossible to carry out in the 
present case due to the limited thickness of the boundary layer. 

A comparison between the measured values and calculated values 
of the phase angle of velocity is shown in Fig. 20 for the 
experiments presented. The two points shown in black which deviate 
from the first bisecting section correspond to the stations where 
at certain instants there is a return flow. Their dispersion 

should not be surprising, since the law of friction loses its meaning 

when the parietal friction reduces to zero. We have also reported 


26 


in this figure the results obtained by M.H. Patel [12] for a 
configuration In which the external flow perturbation is displaced 
considerably by the mean flow velocity: 

Ml « 04 * M/^^flth Q 3r Ifc 


In such a configuration the Instantaneous pressure gradient 
is created both by time derivatives and in X of the external velocity. 
Experimentally, we have considerably larger phase angles which are 
quite predictable by the present calculations. Let us note that in 
the case of the Patel study, there was no decline In the phase near 


the skin. 



Figure 20 - Phase angle 
calculation. Comparison 
with the experimental case. 


4.4. - Integral Method. 

Application to the 
experimental case. 

We have seen that the 
experimental profiles obtained 
at various instants in the cycle 
are correctly represented by a 
family of profiles, with only one 
parameter, deduced from similarity 
solutions. . This permits the use 
of a simple integral method to 
calculate the unsteady flows 
which may lead to separation of the 
boundary layer. Such a method 
has already given good results in 
the case of isolating boundary 
layers on a flat plate configuration 
[3-4]. We shall recall briefly 
its principles. 


r The overall equations of continuity and quantity of motion 
are Integrated step by step over time. These equations are 
written In the incompressible case: 
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it 

1% 




(1*) 


£i 

i 


>• 


u* 9s u> »♦' 


( 15 ) 


6 is the thickness of the boundary layer and Is placed In 
a steady state: 



H* S t /» Cf/t- 


The five unknowns which appear In the system [x#*-*'*-' c f - ft, *. s ) /lh 
are estimated by using additional relationships deduced from the 
similarity solutions already presented. 


CALCULATION- OF. THE DRIVE COEFFICIENT 


This coefficient, having the shape a/i* t appears 

in the overall equation of continuity. It may be calculated from 
the local equation of continuity by making an assumption of similarity 
and by using the small parameter r« {57? which was already 

used for the equation of speed momentum. We have then the 
relationship : 


which permits the expression of vertical velocity V. 


(16) 


At the external zone of the boundary layer we have then: 

s I £>*■ 


( 17 ) 


The drive coefficient Is expressed then from this relation- 
ship and from (5) in the form: 


ii - A. = f * - 

Jx u e V* 


(18) 
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where P Is a function of 0 which results from the similarity 
solutions and which has been represented by the relationship: 


9 • 0,074 « - 1,0*57 / c 


( 19 ) 


In comparison with the steady case. It may be observed that 
an additional term Is added In this drive coefficient : - -1 

The other closure relationships of global equations have been 
obtained in paragraphs 3 .2. Uhls Is the law of friction (10) and 
of relationships (7) and (0) deduced directly from the definition 
of the adverse velocity profiles. It is recalled that these 
relat ionships are : 

»*(•) 

t x y 

G - 2LA 
xr 


The integration of the global equations of momentum and of 
continuity (14) (15) to which the closure relationships may be added 
requires knowledge of the boundary conditions. We must know the 
external velocity as a function of X at each instant and at the 
Initial conditions of the boundary layer: at instant t * 0 it is 

necessary to know the state of the boundary layer as a function of 
X and by giving, for example, the distribution of 6^ and 0. In the 
case of a periodic flow this condition does not affect the result, 
the latter also becomes periodic after a sufficient period of 
calculation. , The other initial condition deals with the given 
boundary layer in X » 0 and at each instant. This assumption is 
important, because in the steady state the accuracy of the estimates may 
depend on it. 


The Integral method which we have described has been applied 
to the experimental case presented. , The external velocity 
la Introduced from the harmonic analysis, limited to the second term, 
of the experimental distribution of velocity. The data of 6 and 8 
as a function of time at the first point of calculation in X are also 
deduced from the harmonic analyslstof integral thicknesses measured. 


Figures 21 and 22 demonstrate the estimated development of the 
displacement thickness 6^ and of the shape parameter H represented 
by their harmonic analysis, , Tie relative amplitudes of the fluctuation 
of H and of 6^ are related to those of the external velocity. 



Key: 


a) Experiment 

b) Unsteady measurement 

c) Quasisteady measurement 

d) Steady measurement of the 
mean flow 
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Figure 21. Harmonic analysis of 

the shape factor. Measurement - experiment. 
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. The measurement does not exceed the station X ■ ^*50 mm. since down- 
stream from this point begins a return flow which Implies the 
transmission of data from the downstream toward the upstream end 
and thus requires knowledge of the additional boundary conditions 
at certain Instants. A more detailed discussion and careful study of 
this problem relating to these boundary conditions Is presented In [5]. 

> Tie evolution of mean values has been determined from a steady 
measurement for the mean flow. The small difference between this 
measurement and the unsteady measurement confirms what has already been 
observed experimentally, namely the small effect of the unsteady 
state on the development of the mean flow, an effect which is due 
essentially to the nonlinear terms of the equations. We must never- 
theless point out that this effect becomes more Important as the 
shape parameter Increases, l.e. when we come near a region with 
return flow. 

> T> determine the Importance of these unsteady effects in the 
configuration under study, we present also a quasi-steady measurement. 

Tils consists of a definition at each Instant of the boundary layer 
development by a steady method but by applying the natural velocity /15 
distribution which exists at this instant. This means that the 
unsteady terms will be omitted in the equations. The developments 
which are thus calculated from the various values deviate consider- 
ably from those measured or calculated by the unsteady 

method. Tils is particularly obvious for the fluctuating values. 

Tie deviations which appear between the unsteady calculation 
and the experimental results relate to the evolution of the fluctuation 
amplitudes or to the phase angle of 6^ and H does not seem very 
large, if account is taken of the difficulty in measuring these values wJ th 
a good accuracy. If it is assumed, for example, that there is an accuracy of 
3? on the measurement of 6^ and 1 % on that of U e , the accuracy obtained 
on the amplitude relative to the variation of 6^ related to that of 
U Q is about 20 %. - lie accuracy of the phase angles of fij and H is 
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Figure 22. Harmonic analysis of 
the thickness of displacement. 
Measurement . experiment . 

Key: a) Experiment 

b) Unsteady measurement 

c) Quasi-steady measurement 

d) Steady measurement of 
the mean flow 


more difficult to determine but 
must be close to 10 or 20 

5 - CONCLUSION 

The experimental study of an 
oscillating boundary layer sub- 
jected to a quite high mean pressure 
gradient resulting In separation 
demonstrates that the development 
of the boundary layer Is not 
profoundly modified by the unsteady 
characteristic of the flow. Ills 
result extends to the new case 
which we have obtained during this 
study of the flat plate configuration, 
always for moderate Strouhal numbers. 
We may nevertheless conclude with 
certainty that this study is not 
completed. In fact, if the measure- 
ments of mean velocity, in the 
sense of overall average, as well as 
of the longitudinal turbulence 
component have been performed up to 
a point where the return flows 
exist and even beyond that due to 
the use of a laser veloclmeter, we 
have available for the present time 
no Information about a value which 
Is as important as the turbulent 
friction <u* v*>. This shall be 
dealt with In the study that follows. 


Comparison of the experimental profiles of velocity with the 
profiles deduced from the similarity solutions shows a good agreement, 
Including for those high values of the shape parameter close to 2. 
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The existence of a logarithmic region is also brought to light. , The 
development by a method of small perturbations of this logarithmic 
law gives Important results relating to the shape of the phase angle 
profiles of velocity in the boundary layer and permits an accurate 
estimate to be made of the maximum phase angle. 

Finally, an unsteady Integral method using the global equations 
of momentum and of continuity together with the closure relationships 
deduced from the similarity solutions, makes it possible to predict 
quite accurately the development of the boundary layer in the con- 
figuration presented. Such a method is nevertheless limited to a 
range where there are no negative velocities, on the one hand because 
the closure relationships are used, and on the other hand because of 
the problem of the boundary conditions. The interest of such a 
calculation method is in its simplicity, which is Important for 
applying practical solutions. It remains now to define its range of 
validity, especially for high values of the Strouhal number. 
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